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The Coulomb drag between parallel, ballistic quantum wires is studied theoreti- 
cally in the presence of a perpendicular magnetic field B. The transresistance Rd 
shows peaks as a function of the Fermi level and splitting energy between the ID 
subbands of the wires. The sharpest peaks appear when the Fermi level crosses the 
subband extrema so that the Fermi momenta are small. Two other kinds of peaks 
appear when either intra- or inier-subband transitions of electrons have maximum 
probability; the m£ra-subband transitions correspond to a small splitting energy. 
Rd depends on the field B in a nonmonotonic fashion: it decreases with B, as a 
result of the suppression of backscattering, and increases sharply when the Fermi 
level approaches the subband bottoms and the suppression is outbalanced by the 
increase of the Coulomb matrix elements and of the density of states. 
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I. INTRODUCTION 



Experimentally jl] and theoretically |2j momentum transfer between spatially separated 
electron layers or Coulomb drag has been studied mostly between two-dimensional (2D) 
layers. Theoretically, this drag has also been studied between very long one-dimensional 
(ID) layers (quantum wires) in which the wire length L is much longer than the mean 
free path (diffussive regime ||) and recently between ID layers of length L <C Zj, in 
which the electron motion along the layer, at low temperatures, is mostly ballistic M || 
(ballistic regime). Even when most of the electrons pass along the wires without collisions, 
a few of them experience backscattering due to the interaction with the electrons of the 
other wire and this modifies the time-averaged distribution functions in such a way that 
the drag effect occurs. In both regimes the drag response is found to be maximal when 
the subbands in the two wires line up precisely. It is important that in the ballistic regime 
the transresistance is determined only by the Coulomb interaction between the electrons 
and such basic properties of the layers as the number of occupied subbands, and does not 
include the relaxation characteristics of the system such as scattering times. Therefore, the 
ballistic regime provides a possibility to obtain more direct information about the Coulomb 
interaction in ID electron systems. 

Motivated by the results of Refs. || and [|] we undertook an extended theoretical study 
of the drag in the ballistic regime, without tunneling between the wires, but in the presence 
of a perpendicular magnetic field B. In Sees. II and III, we generalize the theory of Ref. 4 
to include the effects of intersubband transitions in electron-electron collisions and account 
for the influence of a magnetic field on the Coulomb drag; a limited account of this influence, 
valid when only the lowest subbands in the two wires are occupied, appeared in Ref. [5]. 
Concluding remarks follow in Sec. IV. 
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II. GENERAL FORMALISM 



We use a model of a four-terminal double- quantum- wire system, as shown schematically 
in Fig. 1, similar to the systems investigated in the "directional coupler" problem |J. Two 
closely spaced quantum wires, numbered 1 and 2 and centered at y\ and 1/2, are contacted 
independently to four leads at x = and x = L, where L is the length of the wires. 
The leads have chemical potentials /^^(O) = /i^ 2 , and ^^{L) = /i^ 2 . Applying the bias 
^ = {f^t ~~ \h. )/ e to the leads of wire 2 (drive wire) we obtain the current / flowing through 
it. This current induces a voltage Vd = {f^f — A*r)/ e * n w ^ Te 1 (drag wire). This is the 
typical setup for drag measurements [p]. We assume that the barrier between the wires is 
high enough to allow the neglect of tunneling. 

Below we assume that the electrons in each wire are parabolically confined by the poten- 
tials Uj = e < j+m*fl'j(y — yj) 2 /2, j = 1, 2. In the presence of a perpendicular magnetic field B, 
introduced through the vector potential A = (—By, 0, 0), the normalized wave functions are 

tfinfcfoy) = e* k *x 3 nk(y), Xjnk(y) = (7r 1 /%.2»n!)- 1 / 3! Hntty-Y^/e^eM-iy-yj) 2 /^), 
where n is the ID subband number, k the wave vector of electrons, and H n (x) the Hermite 
polynomials. Neglecting spin splitting the corresponding energy spectrum £j n (k) reads 

Sink = £• + hUj{n + 1/2) + {h 2 /2m 3 ){k - yj /f c f. (1) 

Here to 2 = uj 2 +VL 2 , uj c = eB/m* is the cyclotron frequency, rnj = m*uj 2 /VL 2 , i c = {U/rrfuJc) 1 / 2 
is the magnetic length, £ 2 = %/m*ujj, and Yj = [Vt 2 yj + hu c k/m*]/u 2 are the fc-dependent 
centers of the oscillators. The expressions for the kinetic energies (h 2 /2rrij)(k — yj/£l) 2 of 
the electrons can be simplified by a gauge invariant transformation resulting in a shift of k 
by an arbitrary constant. Since we neglect tunneling, we do not consider electron transitions 
between the wires and can make such shifts independently for each wire; this does not affect 
the kinetic equations written below. Explicitly, we will shift the wave vectors in the manner 
k — y\ji 2 c — > k for wire 1 and k — > k for wire 2. Then the kinetic energies in Eq. (1) 

will read h 2 k 2 /2rrij and the oscillator centers Yj = yj + {huj c /m*uj 2 )k. 
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If the distribution functions fj k (x) = fj k change over distances much longer than both 
the electronic wavelength n/k and the characteristic radius of the interaction potential, we 
can write the Boltzman kinetic equations as 

hk df jnk (x) _ 4tt v v , 2 

— fa ft 1^ 2^ | J % n ;«'nl fc ' K ' 9) °{ £ jnk + Sj'n'k' ~ £ jm,k-q ~ ^j'n'^k'+q) 

1 j'k'q nn'jiiTij 

x [/jnfc(l — fjni,k-q)fj'n'k'(l — fj'n^k'+q) ~ fjm,k-q(l ~ fjnk) fj'n'^k'+qK^ ~ fj'n'k')], (2) 

where the collision integral accounts only for electron-electron scattering. The Coulomb 
matrix elements M^f ^ n , n (k, k', q) are given by 



M nX J n>n( k i k 'i<l) = ~ J d V J d V' K o(\<l\ \V ~ y^)Xjnk(y)Xj'n'k'W)Xj'n' 1 k' +q(y')Xjn 1 k- q (y) , (3) 

where k is the dielectric constant and K the modified Bessel function. 

It is convenient to write separately the distribution functions for the forward- and 
backward-moving electrons as ff^ = f jk \k>o and = f jk \ k<0 , respectively. For 

these functions the boundary conditions are given in the Landauer-Buttiker approach by 
fjnk( Q ) = f( £ jnk - Sfi}) and f~ nk (L) = f(e jnk - Sfij), where Sfif = /if - fi, fi is the equilib- 
rium chemical potential, /(e) = [e( £- ^/ fcsT + and T the temperature. For j — 1 and 
j = 2, Eq. (2) gives two coupled kinetic equations whose solution allows us to express the 
unknown potentials /if and fif through the fixed //J an d lA values and thereby calculate 
the transresistance. 



III. RESULTS 



If most of the electrons move through the wires balistically, Eq. (2) can be solved 
by simple iterations. The zero-order approximation gives fj~ nk (x) = f(sj nk — S/ij~) and 
fjnk( x ) = f( £ jnk — d^j)- Substitution of these functions in the collision integral gives a non- 
zero contribution for backscattering collisions between the electrons of different wires. This is 
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the main contribution which will be considered in detail in the following. If more than a single 
subband in a wire is occupied, the intersubband transitions within one wire also contribute 
to the collision integral of Eq. (2) (if only the lowest subband is occupied, the intralayer part 
of the collision integral completely vanishes because of the relation q = k — k' following from 
the momentum and energy conservation rules). Howewer, within the iterative approach used 
here, we can neglect the influence of the intralayer collisions on the distribution functions of 
the drive layer (j = 2), since the transport regime is nearly ballistic. Further, the intralayer 
collisions do not modify considerably the distribution functions of the drag layer (j = 1) 
because Hi — Hi is assumed to be much smaller than [i\ — iHi an d the main effect on 
fi nk (x) results from the interlayer Coulomb interaction. Considering only contributions 
linear in V, we substitute the equilibrium Fermi-Dirac functions fi nk (x) = f(si nk ) in the 
collision integral and obtain 

fink( x ) = f( e ink - tint) - eVimx/h^Xn^x, (4) 

fink( x ) = /(sin* - Sfii) + eV{mi/hk)\ n (-k)(x - L), (5) 

where the factor 

X n (k) = ^j; 2 { \ M nTn\n'n(k, fc', qtf S(e lnk + E 2 n'k' ~ S lni , k - q ~ E^k'+q) 

B nin'ti'j k'q 

X/(£lnfc)[l - /(£lnifc-g)]/(£2n'fc')[ 1 ~ f (^n'.k'+q)] k , >o k , +q<0 ~ [-]fc'<0,fc'+ 5 >o} ( 6 ) 

is determined by the Coulomb matrix elements and the equilibrium distribution functions 
only. The current flowing in the drag wire is given by 

= - E / dk (^M) ifUx) - A^(x)]. (7) 

Id does not depend on x due to the property J2 n I^°oo ^n(k)dk = 0, which follows from 
detailed balance. Substituting Eqs. (4) and (5) into Eq. (7), using the requirement Id = 0, 
and defining the transresistance Rd as Rd = —Vd/I through the ballistic current / = V/Rn, 
where Rn = h/2e 2 N is the ballistic resistance of wire 2 and iV the number of occupied 
subbands, we finally obtain 



flT, N D~1 N-l o oo q 2 

^ = at- at- 2 , T E E / dk' dq dk\M™J , n (k,k',q) 
NN D e 2 k B T n ni=Qn , K=0 J -co J-k> Jo 

X^(^lnfe + ^2n'k' — £lni,fc-g ~ £ : 2n' 1 ,fc'+g)/(^lnfc) [1 ~ / (^lnifc-g)] f(^2n'k' ) [1 — /(^n^fe'+g)] (8) 

Here N D is the number of occupied subbands in the drag wire (wire 1). Note that the 
introduction of both N and Np assumes that the theory is valid when the Fermi energy 
pi — — hujj(n +1/2) with respect to the highest occupied level is larger than k B T. This, 
of course, implies that the ID subband separations hu 1 and hu 2 are much larger than k B T 
and is true at T ~ 1 K for electrostatically defined electron channels. 

Below we consider the case of identical wires Qi = Q 2 = which entails ui = u 2 = oo, 
h = h = f-i and m\ = m 2 ). To further evaluate expression (8), it is convenient to detach 
the contribution R$ from Rp that expresses the equality n + n' = n\ + n[ for which the 
energy conservation law gives q = k — k' . Then we have R D = R$ + R$ with 



m* k B TLu 



N D -1 N-l 

v E 



5, 



n+n' ,m+n'i 



i 2 , 

n.n. 



NN D h 5 e 2 W n ,tr=o„^ =0 knkn'iK + K>) sinh 2 A n „ ; 



X 



-^ni n| n'ni^n, k n > , fc n + k n ') 



(9) 



Here A nn / = [A + huo(n — n[)]/2k B T and A = ej - e" is the interwire splitting energy 
between the lowest subbands. Further, k n = (u/Q)[2m*(fj, — — hu(n + 1/2)) /h 2 ] 1 ^ 2 and 
fc n / = (w/0)[2m*(/z — £° ~~ huo{n' + i/2))/^ 2 ] 1 / 2 are the Fermi wave numbers for the states 
1, n and 2, n', respectively. The part Rq\ corresponding to n + n' ^ ri\ + n[, is obtained as 



Rd = 2hNN D e 2 k B m 2 J^JfJo dk L ^^-W.W^MM') 



x 6 ( kk! + uj 2 (n + n' - n x - n[) /Q 



2,,2\ 



Ml 2 rf in / n (k, k',q) 



f{^lnk)f{^2 



n'k' 
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x [1 - f((e lnk + £2**0/2 + A(k, k')/2)} [1 - f((e lnk + e wv )/2 - A(k, k')/2)} (10) 

where p(k, k') = [(k + k') 2 /A + u 2 (n + n' - m - n^) /Vt 2 (i 2 } 1 / 2 , q = (k + k')/2 + p(k, k'), and 
A(k, k') = A + hw(ni - n[) - (Q/uj) 2 h 2 p(k,k')(k - k')/m*. The statistical factor in Eq. 
(10) is small unless \£\ n k — fj,\, \£2n'k' ~ ^\) an d \A(k, k')\ are small enough and comparable 
to ksT. This allows the integrals over k and k' are to be carried out in narrow regions 
around k n and k n i, respectively. We used the same property to reduce the contribution Ftp* 
to expression (9). Although the requirement \A(k,k')\ ~ ksT imposes certain restrictions 
on the values of /i, A, and ui, the processes with n + n' ^ m + n[ can give a considerable 
contribution to Rd, especially for \ A + Tioj(n — n' 1 )\ ^> ksT and small. We stress that the 
previous theoretical work Q on the Coulomb drag in the ballistic regime took into account 
only the processes with n = n\ and n' = n[, thus neglecting other processes completely from 
the beginning. The numerical calculations given below demonstrate that this limitation is 
considerable in many cases. 

If only the lowest subbands are occupied in each wire, i.e., for n — n\ — n' — — 0, the 
calculation of the transresistance is considerably simplified. Only Rff contributes to Rd 
and Eq. (9) can be rewritten as 



Expression (11) is convenient for assessing the magnetic-field dependence of the transre- 
sistance Rd ■ It directly demonstrates a significant reduction [5] of the drag effect by 
the magnetic field B, mostly due to the exponential factor. The decrease of Rd starts as 
Rd{B) — Rd{0) ~ —B 2 and becomes exponential with increasing B. The physical reason 
for this decrease is the suppression of backscattering in electron-electron collisions as the 
oscillator centers for forward- and backward-moving electrons are pulled apart by the mag- 
netic field. The characteristic field B for this suppression depends on the position of the 



Rd 



2e 2 m* 3 uj 6 Lk B T (A/2k B T) 2 



nf^K^hhrih + k v ) sinh 2 (A/2A; B r) 




(11) 
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Fermi level e F and is estimated as B ~ (m* / 'e)(m*£l 3 /Ti) 1 / 2 / '(ki + ky). If £p is not far 
from the subband bottoms, B is big and the suppression is weak. If e F is well in between 
the ID subbands, B is estimated as 1 tesla for typical wire parameters. However, when 
Ef, with the increase of B, approaches the subband bottom, the opposite effect takes place: 
the transresistance increases because the wave vectors hi, ky, and q become progressively 
smaller and the suppression of backscattering becomes less important than the increase of 
the Coulomb matrix element and of the density of states. 

Below we present numerical results for the transresistance Rp, expressed in units of the 
fundamental resistance Ro = h/2e 2 , at T = 1.3 K, L = 0.4 /mi, d — \yi — t/2 1 = 50 nm, 
h£l =4 meV, m* = 0.067m , and k = 13. Figure 2 shows the dependence of R D on the 
Fermi energy e F) defined as Ep = A* — {e\-\-e^)/2, calculated at B = and B = 1 tesla. The 
calculations were done assuming that up to two subbands can be populated in each wire. 
As seen in part (a), for A = there are pronounced sharp peaks of Rd when e F crosses 
the bottoms of the first and second subbands. The sharpness of the peaks is explained by 
a strong enhancement of the Coulomb collision probability when k n , k n i, and q are small, 
cf. Eqs. (9) and (10). In Fig. 2 (b), for A = 1 meV, one can see three peaks; the middle 
one appears after ef crosses the bottom of the second subband in wire 2. This peak exists 
due to the processes with n + n' ^ n\ + n' v The third, most prominent peak in Fig. 2 (b), 
appears after ef crosses the bottom of the second subband in wire 1, so that two subbands 
in both wires are populated. The processes with n + n' ^ n\ + n[ give the main contribution 
to this peak as well. 

The application of the magnetic field shifts the peaks to higher Fermi energies, due to the 
increased confinement energy, and sharpens them due to the suppression of backscattering 
in the regions far from the subband edges. The resulting decrease of the transresistance 
due to this suppression is illustrated in Fig. 3 for one (ef = 3.5 meV) and two (ef = 8 
meV) subbands populated. These dependences are non-monotonic: when, with the increase 
of B, £p approaches the first or second subband bottom, Rr> starts to grow sharply. This 
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dependence is basically the same for both A = and A = 1 meV. 

Figure 4 shows the dependence of R D on the level splitting A at several constant values 
of /i — el = Sf — A/2. This means that the subband positions of wire 1 remain constant with 
respect to the Fermi level but those of wire 2 do not; this can be experimentally achieved, for 
example, by changing the voltage of the gate adjacent to wire 2 while keeping that adjacent 
to wire 1 at a constant voltage. The curves are plotted for one (a) or two (b) subbands 
populated in wire 1 but for different /x, far from (solid) and close to (dashed) the upper 
populated subband edge. Both curves of Fig. 4 (a) show two peaks: the sharp ones appear 
when the second subband of wire 2 becomes populated while the broad ones appear when 
the second subband of wire 2 is aligned with the first one of wire 1, at A = hu ~ 4.36 
meV. A similar behavior is seen in Fig. 4 (b). At large negative A only one subband is 
populated in wire 2 while at A ~ —1.5 (solid) and —0.5 meV (dashed) the second subband 
of wire 2 becomes populated as well. This transition is reflected by strong and sharp peaks 
in R D . Other strong peaks appear at A = 0, when the subbands are aligned; note that 
on the dashed curve such a peak merges with that at A ~ —0.5 meV and is not resolved. 
The minor peaks in the regions of negative and positive A exist due to the inter-subband 
transitions with {n, ni, n', n^) = (0, 1, 0, 0) and (1, 0, 0, 0) and (n, ni,n', n' x ) = (0, 1, 1, 1) and 
(1,0,1,1), respectively. The maxima of these peaks occur when A(k,k') ks A(k n ,k n >), cf. 
Eq. (10), goes to zero. Thus, the level-splitting dependence of R D shows a rich structure of 
peaks indicating that both the intra- and mter-subband transitions of electrons contribute 
to Rd- 

All calculations described in this section were repeated for different values of the interwire 
separation d. An increase of d considerably decreases the transresistance: Rd drops by more 
than one order of magnitude as d varies from 40 to 60 nm, mainly due to the dependence of 
the Bessel function on its argument. However, all qualitative features presented above are 
preserved. 
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IV. REMARKS AND CONCLUSIONS 



The treatment of the drag effect in the ballistic transport regime demonstrates the salient 
properties of electron-electron collisions in double-layer quasi-lD electron systems. The re- 
duced dimensionality dramatically decreases the scattering probabilities at low temperatures 
due to the restrictions imposed by the momentum and energy conservation laws. As a re- 
sult, the transresistance shows peaks as a function of either the Fermi level position or the 
interlayer level splitting energy. The peaks always appear when the Fermi level crosses the 
bottom of a subband, so that a new subband n is involved in the scattering prosess; the 
Fermi wave number k n for this subband is small, the density of states is high, and this results 
in a higher scattering probability. When subband n is aligned to another one, the conser- 
vation rules allow electron transitions inside the subband n, the corresponding momentum 
transfer hq ~ 2hk n is small, and the Coulomb matrix element is large thus giving rise to an 
additional increase of the peak. Next, the peaks appear when two subbands from different 
layers are aligned; this favors transitions which conserve the sum of the subband numbers, 
n + n' = ni+n'i, especially the transitions between the electrons inside the aligned subbands, 
cf. Eq. (9). Finally, peaks occur under special conditions, for A(k n ,k n ') ~ 0, cf. Eq. (10); 
this implies a maximum probability for intersubband transitions with n + n' ^ n± + n[. 
Although the peaks associated with these transitions are usually weaker than those under 
subband alignment, they give a considerable contribution which cannot be neglected. The 
described rich structure of the peaks is best seen in the level-splitting dependence of the 
transresistance shown in Fig. 4. 

A magnetic field B applied perpendicular to the wire plane reduces the overlap be- 
tween the wave functions for forward- and backward-moving electrons and thereby tends to 
suppress electron-electron scattering. This results in a decrease of the transresistance. In 
addition, the application of B modifies the quantization energies and leads to a shift of the 
subbands with respect to the Fermi level. Since the scattering probability increases when a 
subband edge comes close to the Fermi level, the transresistance Rp may increase with the 
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increase of B. Therefore, the dependence of Rd on B is basically non-monotonic as shown 
in Fig. 3. 

The results obtained here are valid when the ID electron gas in either wire is described 
as a normal Fermi liquid. We used this model because the wires are short, the transport 
is nearly basllistic, and the properties of the ID electrons are determined by those of the 
2D reservoirs they are injected from. The case of the Coulomb drag between ID electron 
systems described as Luttinger liquids has been studied in Ref. 7. 

Concerning experimental results we are aware only of those of Ref. || where the tran- 
sresistance Rd was measured as a function of side gate voltages controlling the confining 
potentials of the parallel, submicron-long quantum wires, thus allowing to change the posi- 
tions of the ID subbands with respect to the Fermi level, the interlayer subband splitting A, 
the wire widths Wj, and the interwire distance d. It was found that Rd shows sharp peaks 
when the Fermi level crosses the bottom of a ID subband. When the gate adjacent to the 
drag wire was kept at a constant voltage, correponding to one populated subband in it, the 
transresistance, as a function of the voltage of the gate adjacent to the drive wire, showed 
two peaks. The shape and position of these peaks permit us to identify them with those of 
Fig. 4 (a), since the situation described by Fig. 4 (a) corresponds roughly to this type of 
measurements. These experimental results provide qualitative support for our theoretical 
predictions. However, as no formal connection is made in our model between the gate volt- 
ages and the parameters ep, A, d, and Wj, we cannot attempt a more detailed comparison. 
Such a connection requires a detailed knowledge of the gate-induced modification of the 
double-wire confining potential which could be obtained only by a self-consistent solution 
of the electrostatic problem for the three-gate structure investigated in Ref. 8. We expect 
though that further experimental and theoretical work will test sufficiently the drag in the 
ballistic regime and our results. 
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FIGURES 



FIG. 1. Schematic diagram of a coupled quantum-wire device. 

FIG. 2. Dependence of the transresistance Rd on the position of the Fermi level for (a) aligned, 
A = 0, and (b) shifted, A = 1 meV, levels in quantum wires. The dashed and solid curves 
correspond to B = and B = 1 tesla, respectively. The other parameters are listed in the text. 

FIG. 3. Dependence of Rd on the magnetic field B, with one (ep = 3.5 meV) and two {ep = 8 
meV) populated ID subbands, at A = (a) and A = 1 meV (b). 

FIG. 4. Dependence of Rd on the level splitting energy A at B = 1 tesla when one (a) or two 
(b) populated subbands of wire 1 remain constant with respect to the Fermi level, e^ — A/2 = const. 
(a): e F -A/2 =3.5 meV (dashed) and 4 meV (solid), (b): e F -A/2=7 meV (dashed) and 8 meV 
(solid). 
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Fig. 1 . O. E. Raichev et al. , Coulomb drag. 
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Fig. 4. O. E. Raichev etal. Coulomb drag. 



